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L~:> , We study the mean traversal time r for a class of random walks on Newman- Watts small-world 

' networks, in which steps around the edge of the network occur with a transition rate F that is 

^jq' different from the rate / for steps across small-world connections. When f ^ F, the mean time r 

, to traverse the network exhibits a transition associated with percolation of the random graph (i.e., 

small-world) part of the network, and a collapse of the data onto a universal curve. This transition 
was not observed in earlier studies in which equal transition rates were assumed for all allowed steps. 
QQ ] We develop a simple self-consistent effective medium theory and show that it gives a quantitatively 

correct description of the traversal time in all parameter regimes except the immediate neighborhood 
of the transition, as is characteristic of most effective medium theories. 
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INTRODUCTION 



There has been intense recent interest in statistical and dynamical properties associated with random networks 
possessing so-called "small world" (SW) properties, i.e., networks whose bonds possess a large degree of local clustering, 
but with a relatively small minimal path length connecting nodes of the system. At least some of this interest has 
arisen from the realization that algorithmically constructed small-world networks (SWNs) introduced by Watts and 
Strogatz 0], and by Newman and Watts appear to share many statistical features with social networks, and thus 
^ . form a useful topological substrate on which to model dynamical processes relevant to the medical and social sciences 
O I [3- Power grids and information networks may also possess properties associated with small- world or scale- free 
networks 3, and even certain aspects of polymers may be understood in terms of an underlying small- world structure 
i- 

Independently of their possible applications to a wide range of problems in the social, physical, and information 
sciences, small-world networks and their extensions are inherently interesting because they provide a class of disordered 
structures similar to those that have been studied in the condensed matter literature for many years |E S H S ^| ■ 
' There is an extensive literature, for example, in which the dynamics of random walkers diffusing through various kinds 
00 1 of disordered media have been studied throurii a var iety of methods, including numerical simulations 0, percolation 
O ; ideas 0, effective medium theories (EMTs),]lH 0,113 and others. 

1^ . In the present paper we study traversal times associated with random walks on the Newman- Watts small-world 
' networks (NW-SWN's). These structures are defined on a ring of N sites, each of which has bonds that connect it 
. to its 2K nearest neighbors. On this translationally-invariant lattice a (generally) disordered small-world structure 
' is imposed, by associating with probability q a SW connection along each of the remaining N{N — 2K — l)/2 bonds 
I , in the system. For a given value of g, there are on average Ugw = q{N — 2K — l)/2 qN/2 small-world bonds per 
■ site in the system, a parameter that allows for a useful comparison between SWNs of different size. Except for the 
^ r endpoints at g = 1 and q = 0, where the system is perfectly ordered, the network forms a disordered system on which 
various forms of dynamics can be studied. 

The new feature of our analysis, which is based upon a master equation for the site occupation probabilities P„ (t), 
is that we take the jump rate / associated with SW connections to be independent of the jump rate F associated with 
steps around the edge of the network. Models of this sort might arise, for example, in attempts to design, modify, 
' or optimize existing information networks in order to reduce the mean access time, by incorporating a small number 
of fast connections into an existing random network already possessing a large number of other (perhaps slower) 
connections. By introducing this simple extension to the situation in which all allowed steps occur with the same 
probability (and hence same transition rate), our theory is able to reveal interesting universal behavior that occurs 
for f It including a percolation transition that occurs at a critical fraction of SW bonds. Such a transition, which 
is clearly a consequence of a well-known result by Renyi and Erdos ^3 for random graphs, was not observed in earlier 
random walk studies of the NW-SWN [l^ll^ . in which only equal jump probabilities (/ = F) were considered. The 
transition that we observe numerically and analyze using a simple self-consistent effective medium theory, is critically 
dependent on the fact that the topological disorder associated with the small- world connections is "quenched" . We 
show, for example, that an exactly solvable model with "annealed disorder" , in which the walker decides at each step 
whether to move along the edge or to a randomly chosen site elsewhere in the network, shows no such transition as a 
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function of the branching ratio between steps around the edge of the ring and steps associated with smah- world type 
"shortcuts" . 

The rest of the paper is laid out as follows. In the next section we introduce a master equation description of 
continuous time random walks on Newman- Watts small-world networks, and review the relation between the Green's 
functions associated with those equations to the mean time r for a random walker on such a network to traverse 
the system. We then present numerical calculations of the mean traversal time and show that for f ^ F there is a 
transition in the underlying transport mode, between motion that is dominated (and limited ) by edge diffusion and 
motion that takes place predominantly among small-world connections in the system. In that section also we consider 
a model with annealed disorder and show that it does not display the transition exhibited by the model with quenched 
disorder. In Sec. ^we analytically consider the traversal time for NW-SWN's using a simple self-consistent effective 
medium (EMT) theory similar to that which has been used to understand other disordered systems. As in previous 
work, we find that EMT provides an accurate and computationally efficient means for calculating transport properties 
as a function of the underlying parameters of the system, except in the immediate vicinity of the transition where the 
underlying transport mechanism changes abruptly. We are thus able to use it to reliably investigate, in regimes away 
from such points, the manner with which various transport properties scale with system size and connectivity. 



MASTER EQUATIONS, TRAVERSAL TIMES, AND PERCOLATION 

We consider a continuous time random walk with symmetric, translationally-invariant jump rates Fn 
governing transitions between sites n and to on a 1-dimensional ring of N sites. Although it is not crucial in the 
analysis that follows, for specificity we assume that Fn^m — F is constant and non-zero only between a given site 
and its 2K closest neighbors on the ring, which we will refer to collectively as "neighbor sites" . The small-world 
structure of the system is imposed by associating with probability q a nonzero SW hopping rate / along each of the 
N{N — 2K — l)/2 non-neighbor bonds in the system. For q = I and q — 0, the system is perfectly ordered; in the 
former, all sites are connected, in the latter, connections occur only between sites on the periphery of the ring. For 
intermediate values of g a particle moving along the ring encounters sites from which it can move across the system 
without taking steps to its 2K nearest neighbors. The small- world connections thus superimpose a "random graph" 
structure on the otherwise translationally invariant lattice. 

Random walks on any single NW-SWN constructed in this manner can be described by the master equation 

dP,n ^ ^ 

^ ^ F (^Pm-\-n Fm) ^" ^ ^ fmn (-^m+n -^m) (1) 

for the probability Pm(t) for the particle to be at site to at time t, where we assume periodic boundary conditions so 
that Pm+N = Pm- In the first sum on the right describe steps around the edge of the ring, while the last term 
describes SW transitions. In this last term, the f„m are symmetric random rates that in any given realization equal 
/ with probability q and equal zero with probability 1 — g. 

A full solution to the problem involves obtaining the set of propagators, or Green's functions gm.n{t), the solutions 
to Eq. for each initial site n at which the particle may start. These quantities, or their Laplace transforms 
gm,n (s) = /o ffm,n {t) e~'^dt, are readily computed for moderately sized systems N < 10'^ using standard numerical 
techniques. For example, rewriting the equations of motion (^) in the form 

^ + 5]A„„P„(i) = (2) 

n 

which along with Eq. implicitly defines a transition matrix A with elements Amn-, we can write the Green's 
function in the time domain as the m-n element of an exponential 

5™.nW=[e-^*]„,„ (3) 
of the matrix A. Similarly, we can express its Laplace transform as the corresponding element 

5m,n(£) = [(f + A)-'l (4) 

of the resolvent matrix {e + A) ^, which is easily computed using matrix inversion routines. 





N= lO' 



l/N 



10-' N 




W'N Nil 



Mean Number h^,,, of SW Connections Per Site 



FIG. 1: Mean traversal time r as a function of the mean number Usw of small-world connections per site on a SWN of lO'^ 
sites, with values of f /F as indicated. 

Numerical solutions of this sort can be used to find information about the evolution of the probabilities Pm {t) , 
their moments, or other quantities that characterize transport in the system. In the condensed matter literature, 
e.g., much attention has been placed on calculating the diffusion constant D which in an infinite Euclidean network 
characterizes the asymptotic linear growth of the mean square displacement {n^ (t)) ^ Dt of an ensemble of random 
walkers. Insofar as we are interested in the properties of finite networks, in which the mean-square displacement 
always saturates at long times, we focus here on calculating properties that may be more relevant to applications for 
which small-world concepts are currently employed. Specifically, we focus here on calculating the mean traversal time 
r, which we define to be the earliest time, on average, that a random walker visits the point on the ring the farthest 
from where it started. For a walker starting at site n at t = 0, the mean time t,„_„ to arrive at an arbitrary site m 
for the first time is the first moment of the probability density Fm,n{t) for a walker to first arrive at site m at time t 
for these initial conditions, i.e.. 



where in the second line we have expressed the result in terms of the Laplace transform F,n.n (&) of Fm,n (t)- The 
first passage probability density F^.n {t) is related to elements of the Green's functions described above through the 
relation 



which physically expresses the probability for the walker to be found at m at the current time t, in terms of its 
probability to have arrived at that site for the first time at some earlier moment i', and the probability that it is now 
at that site given that it did so. From the convolution theorem for Laplace transforms it then follows that 




(5) 



lim 



dF„i^n (e) 
de 



(6) 




(7) 




(8) 



so that we can write Eq. Q as 




(9) 



The mean first passage time is thus readily computed from a numerical solution to the Laplace transformed Green's 
function of the system Q, or equivalently, through the resolvent of the matrix A. 
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In Fig. ^ we display the results of a numerical calculation of the mean traversal time 

r = iV-i^(r„.+ «,J (10) 

m 

scaled by the average nearest neighbor hopping time F~^, averaged over an ensemble of 100 K — 1 NW-SWN's with 
A'' = 10^ as a function of the average number Usw of small-world bonds per site, for different values of the ratio 
f /F. As expected, the mean traversal time t decreases monotonically as the number Ugw of small- world connections 
increases. In the limit in which Usw < N~^, i.e., in the absence of any small-world shortcuts, the mean traversal time 
reduces to the mean time tq that it takes for the walker to diffuse around the ring, for which 2Fto ~ (7V/2)^ . Also, 
as we might expect, for very small values of f / F, only modest decreases in the traversal time occur with the addition 
of small world network connections. 

For values of / > F, however, the situation appears very different, there generally being a very strong decrease in 
the traversal time near Ugw 1 — 5, and a strong collapse of the numerical data onto a single curve for values of 
Usw < 1- In this region SW connections are sparse, but fast. For Ugw < 1 and f ^ F, the SW connections act as 
short-circuits, but transport across the system is still limited by diffusion within segments of the ring that are free 
of shortcuts, because there is not in this regime a percolating path of small- world connections spanning the system. 
However, as shown by Renyi and Erdos [lll |. in a completely random graph a percolating path develops (as N — > oo) 
as Usw approaches one from below. Thus, as f /F oo, and A'' — > oo, percolation of the random-graph part of the 
NW-SWN network leads to a critical change in the traversal time near ~ 1- 

It is interesting to note that other transport properties of the NW-SWN's, including scaling properties of the 
diagonal element gnn (e) of the average green's function|0|, and the distinct number of sites visited S (t) by a walker 
on the networkjl^ have been studied previously, but only for the case in which f = F, i.e., in which hops associated 
with small-world connections occur with the same transition rate as steps along the periphery. As we see from the 
current study, however, the case F = f lies at the edge of a parameter regime in which the behavior of the system, as 
a function of the number of small- world connections, changes drastically. 

It is natural to ask whether the quenched disorder occurring in the small-world networks is an essential element for 
the observed transport threshold to occur, or whether a simpler annealed disorder model with uniform trans-network 
connections could exhibit a similar transport transition as a function of their strength. We show explicitly below that 
a simple annealed disorder model does not exhibit such a transition. Specifically, we study an exactly soluble model 
in which a walker at any site takes steps to its neighbors on the edge of the ring with rate F, as before, and takes steps 
to one of n randomly chosen non-neighbor sites on the ring with rate /, but in such a way that the n randomly chosen 
sites to which it may move are not fixed in time, but are selected anew each time the particle visits the site. Thus, in 
this model, if we take n = 2nsw the branching ratio between steps associated with short cuts across the system and 
those around the edge is, on average, the same as in the NW-SWN model with quenched disorder already considered. 
Since the a priori probability of the particle moving to any other site on the ring other than its neighbors is the same, 
the dynamics on this annealed disorder model is equivalent to one in which the short cuts allow the walker to move 
to any one of its non-neighbors with a uniform rate 

Lo = 2ns^f/N„, (11) 

where Nn — N ~ 2K — 1 is the number of non-neighbor sites to which it can move. Thus in any given site, the 
walker moves to one of its neighbors with probability 2KF/{2KF + 2nswf) and to a non-neighbor with probability 
Uswf / {2KF + 2nswf)- The resulting master equation for this model, 



F{Prn+n-Prn)+ ^ iU {P,n+n - Prn) (12) 



dt 

is translationally invariant and readily solved by introducing Fourier transformed probabilities P'^ (t) = J^n Pn (t) e**^". 
We find, e.g., that for this annealed disorder model the Laplace transformed Green's functions, g:^ (e) = g:^ Q (e) , are 
given by the relation 

^-(^)- 7Ve(f+iVc.) +^-(^ + ^-) (13) 

in which 



Mean Number n of Short Cuts Per Site 



FIG. 2: Mean traversal time r as a function of the mean number n of trans-network connections per site for the annnealed 
disorder model on a system with 10* sites, with values of f /F as indicated, showing the failure of the annealed disorder to 
capture the transport transition observed to occur on real small-world networks. 



K 

Afc = ^2(1 -cosfcn), (15) 

n=l 

and the sum in p4|l is over all wavevectors k = 27r?T,/iV in the Brillouin zone. 

In Fig. 121 we show how the mean traversal time r scales in the annealed disorder model as function of n (now 
considered a continuous parameter), for various values of / for a system with N ~ 10'*. Clearly, although the limiting 
diffusive behavior seen in the quenched disorder model of Fig. ^ as ^ and / — > is the same as that observed 
in Fig. 12 the previously observed behavior for f ^ F, and in particular the transition near n^u; 1 is completely 
lacking in the model with annealed disorder. Thus, as we might expect, the underlying percolative behavior is clearly 
a property associated with a system with quenched disorder. We note in passing that the annealed disorder model 
described above can be viewed as an approximation to the quenched disorder model in which the random matrix A 
is replaced by its ensemble average (A), i.e., in which 

((e + A)-i)«(£ + (A})-\ (16) 

For disordered systems, approximations of this type are known to give results that are often qualitatively very bad, 
particularly in percolative systems, where they typically wash out any transition that occurs. It is not surprising, 
therefore, that the two models have very different behavior. 

It is clear that the annealed disorder model fails to capture the essential features of the transport transition that 
occurs in the NW small world networks. However, the idea of replacing a disordered system with a translationally- 
invariant one that captures, on large time and length scales, the macroscopic transport properties of the ensemble of 
systems that it replaces is venerable and has a long history in the condensed matter literature. Indeed, the goal behind 
a significant body of theoretical work on various kinds of disordered systems has been to construct an appropriate, 
self-consistently determined effective medium, whose properties captures qualitative features of the actual system 
under consideration. In the next section, we explore this idea, and introduce a simple effective medium theory that 
does, in fact, provide an excellent quantitative prediction of the traversal time for the NW-SWN in all parameter 
regimes except the immediate neighborhood of the critical point. 

SELF-CONSISTENT EFFECTIVE MEDIUM THEORY OF THE TRAVERSAL TIME 

Theoretical justification for the search for a translationally-invariant effective medium whose properties capture the 
essential features of the ensemble of disordered systems that they replace (and which capture the large time and length 
scale properties, typically, for any member of the ensemble) lies with the fact that average transport properties of the 
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FIG. 3: Effective medium parameter wo as a function of tlie mean number Uaw of small-world connections per site on a SWN 
of 10* sites, with values oi f / F as indicated. Inset: Effective medium parameter as a function of the number A'^ of sites in 
the network, for the case f = F. 



ensemble are, in fact, translationally invariant. For example, the ensemble-averaged probabilities pm (t) — {Pm{t)) 
associated with any fixed initial condition evolve in a translationally invariant way, and therefore obey translationally 
invariant equations of motion which, if we knew what they were, would serve to define the effective medium that 
we seek. Indeed, the goal of effective medium theory is to self-consistently determine properties associated with 
this average evolution. To this end, we consider the simplest set of linear, homogenous, and translationally-invariant 
equations of motion that reflect the structure of the original equations, and the symmetry properties of the underlying 
ensemble. In particular, the Laplace transforms Pm (e) of the average probabilities Pm (t) we take to obey the equations 

K K 

ePm (e) - Pr,i (0) = ^ {Pr,i+n - P,n) + ^ W (e) {pr,i+n - Pm) ■ (17) 
n=-K n^-K 

Here, w (e) is a frequency-dependent rate, equivalently, a memory function in the Laplace domain |l4|. connecting 
pairs of sites on the network capable of being connected, in any realization, by a small-world rate /. In these effective 
medium equations of motion, transport around the ring edge is characterized by the same rate F that obtains 
throughout the ensemble, but w (e) must be determined from self-consistent considerations. Note that these EMT 
equations of motion have exactly the same form as the Laplace transform of those describing the annealed disorder 
model of the last section, except that in that model the rate w (e) = to was simply set equal to its ensemble average, 
while in the current treatment w (e) is to be determined self-consistently through other considerations. 

Thus, the theoretical tasks are two: (i) the self-consistent determination of w (e) as a function of N and q (or risw), 
and (ii) the determination of transport properties arising through the solutions to (|17l) . Both tasks require calculation 
of the effective medium propagators gf^ (e), i.e., the solutions to ifTTIl for a walker initially at the origin, pm (0) — Sm,o- 
Since the structure of the equations of motion (|17|l are the same as those of the annealed disorder model, the form of 
the effective medium Green's function associated with ((T7j) is the same as ((T^ , with lo replaced by w (e) : 

Nw 

where 

Using (fl^ we now develop a self-consistent expression for w{e). We proceed as in other EMTs [ElESEl) by 
embedding in the effective medium a typical "fluctuation from the average" associated with the disorder, and then 
requiring that upon averaging over the distribution associated with this fluctuation we recover the effective medium 
in which it was embedded. Thus, we consider the bond between the origin and some other non-neighbor site n on the 
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FIG. 4: Average time to first traverse the network as a function of the mean number Usw of small world connections per site, 
on a SWN of 10'^ sites, with values of f/F as indicated. Solid lines are from EMT, data points are the same as in Fig.0 



ring, and replace the effective medium bond w (e) between those two sites with an actual rate /o,„ = /n (= or /), 
drawn from the ensemble. For a particle placed at the origin, the resulting pmi^) obey 

K K 

SPmis) - S,n,0 ^ ^ F{pn+Tn~Prn)+ ^ W (s) {pm+n - Pni) - ^ni (s) (20) 
n=-K n^-K 

where Am (e) — (<5m,o — Sm,n) [fn — w) {po — Pn)- The solution to the set of equations (j2()|l can be written in terms of 
(Pt as 

1 + 2a„7„ 

where a„ (e) = /„ - w (e) , and 7„ (e) = (e) - {e) . 

We now impose self-consistency, and require that, upon averaging (|21() over the dichotomous distribution 

p{fn) = {l-q)S{fn)+q5{fn-f) (22) 

of the rates /„, on the one hand, and the location n of the small- world connection, on the other, we recover the 
propagator for the effective medium. For this to occur, the second term on the right hand side of H21|l must average 
out to zero. This leads to the self-consistent condition 



N/2 2 /I ^ ~~2 

q[f -w)-/n sr^ (1-9) i«7„ 



^ l-|-2(/-w)7„ ^ l-2w7„ 

n=K+l ' n=K^\ 



(23) 



from which w (e) can be determined. For q = 1 and g = 0, respectively, the solutions reduce to the exact results. For 
g ^ or 1, Ea. (|23|l can be solved numerically without having to perform large matrix inverses. 

Because we are primarily interested in properties of the network associated with large length scales, and thus with 
times much longer than that required for a single hop, we are interested in the Laplace domain on the behavior of 
the system for small values of the Laplace variable e as in Eq. 10. Accordingly, we present in Fig. Ola plot of the 
zero-frequency effective medium parameter wq — lime^o w (e) as a function of the number Hsw — qN/2 of small-world 
connections per site, for different values of the SW hopping rate /, on a NW-SWN with N = 10** sites. On the 
right hand side of the figure, at sufficiently large values of Usw for any value of /, the effective medium parameter wq 
appears well-described by the formula wq qf 2nswf /N. As we move to values of Usw ~ 1 and less, there is a 
convergence of the curves with f > F, similar to that which appeared in Fig.^for the traversal time of the NW-SWN, 
thus giving hope that the self-consistent treatment introduced here provides a qualitatively correct description of 
transport in this disordered system. After some analysis of the results, we find that in the region 1 > nsw^ where the 
different curves with f > F converge, the parameter wq is independent of / and well-characterized by the relation 
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FIG. 5: Mean traversal time as a function of the number A'' of sites in the network, for f — F. Dashed hue is the hmiting case 
for pure diffusion, sohd hues are the results of EMT. Dotted lines indicate scaling relationships discussed in the text. 



wq ~ Nq^F ~ in^^F/N, depicted as a dashed line in the figure. As iV oo, the region of the main part of Fig. 
13 associated with this universal behavior extends downward and to the left, so that for any / there will be a value 
n-swif) below which wq will scale as An1^F/N . In the inset of Figure 01 we show how this scaling relation is obeyed 
as a function of system size for large N and small n^^j, for the case in which f — F. Dotted lines in the inset indicate 
the relation wq ~ An^^F/N, which the curves with small Usw approach for large N. At moderate values of Ugw > 1, 
scaling crosses over to the form wq ~ 2nswf /N, represented by the dot-dashed line, shown for Ugw = 2.5. For very 
small values of ns^, < N~'^, uninteresting deviations from SW scaling occur in the regime where most members of 
the ensemble have no SW connections. 

This interesting behavior of the effective medium parameter wq leads to two important questions. Does the behavior 
seen in the parameter wq lead to a corresponding effect on actual transport observables? And do the predictions of 
EMT for those observables describe the actual system of interest? At least in terms of the mean traversal time the 
answer to both of these questions is affirmative. The derivatives and small-e limit in Q can be explicitly taken for 
the effective medium propagators HI8|I . leading to the simplification 

t™ - N [Go {Nwo) - G,„ {Nwo)] • (24) 

In Fig. 01 we plot the mean traversal time r = rjv/2 as a function of Usw for a system with N = 10'^ sites. Solid 
lines in this figure are predictions of EMT, while the data points are the same as those appearing in Fig. ^ In Fig. 
01 we see that the effective medium theory predicts a transition in the vicinity of Usw = 1 , as well as a collapse of the 
numerical data for large f /F onto a single curve for Ugw < 1- For all values of /, traversal times saturate for very 
small to a value such that 2Ft (iV/2)^, consistent with pure diffusion around the edge of the ring. In the region 
1 ^ nsw ^ where the data collapse occurs for f > F, the traversal time of the numerical data and the EFT are 
both described by the functional relation r ~ {2qF)^^ = N/Aug^F, which is a reflection of the mean number of steps 
Ft ^ {2q)^^ the particle must take along the edge of the ring before it encounters a SW shortcut across the system. 

Comparison of the data and the solid curves in Fig. ^ shows that, as in other EMTs, the theory derived above 
appears to be numerically accurate for all parameters except those near the underlying percolation transition, the 
effect of which only becomes apparent for large values of //F in the neighborhood of ~ 1- Provided that we avoid 
this particular regime, i.e., the region bounded by the dashed curves in Fig.^ we can reliably use our EMT to calculate 
accurately other properties of the system. The two dashed curves in that figure represent limiting cases for very large 
f/F> 10® for the EMT (on the left) and the numerical results (on the right). The discrepancy between EMT, and 
the exact numerical results does not decrease significantly with increasing TV. Thus, the EMT of the present paper 
is typical of those for other disordered systems, which often capture the essential behavior of the system, but fail to 
accurately reproduce critical properties. Of course the main advantage of the EMT we have derived is that we can 
now extend its calculations to larger networks (i.e., larger values of N) for which numerical solutions to the individual 
equations of motion become prohibitively difficult. As an example, we present in Fig. [3 a plot of traversal time t 
as a function of system size N. This plot has been computed using EMT for fixed values of Usw, for the case f — F, 
which is accurately described by EMT for all rigw, as is evident from Fig. 01 Clearly, the predicted scaling r ~ iV, 
shown as dotted lines in Fig. jSl obtains for any value of Hstv for sufficiently large values of A^. Indeed, for small Usw, 
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the time to traverse the system fohows closely the diffusive result t ^ N"^ until the system size is sufScicntly large 
that there is a significant probability to encounter a shortcut before reaching the other side of the ring via diffusion 
around the periphery. 

SUMMARY 

In this paper we have considered the mean traversal time r for a class of random walks on Newman- Watts small- 
world networks, in which steps around the edge of the network occur with a transition rate F that is different than 
the rate / for steps across small-world connections. Using numerical calculations of the Green's function we obtained 
present numerical data for the mean traversal time and show that for f ^ F there is a transition in the underlying 
transport mode, between diffusion limited and connection limited motion. We showed that a model with annealed 
disorder does not display a transition reflective of the underlying percolation transition occurring on the random 
graph part of the network. We then developed a self-consistent effective medium (EMT) theory for the traversal time 
that provides an accurate and computationally efficient means for calculating this quantity except in the immediate 
vicinity of the transition. As a result, we have been able to use EMT to determine how various transport properties 
scale with system size and connectivity on the Newman- Watts small world networks. 

This work was supported in part by the National Science Foundation under Grants DMR-0097204, DMR-0097210, 
INT-0336343, and DARPA-N00014-03- 1-0900. 
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